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Thermal Analysis on Plume Heating of the Main Engine on the 
Crew Exploration Vehicle Service Module 

Xiao- Yen J. Wang and James R.Yuko 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

The crew exploration vehicle (CEV) service module (SM) main engine plume heating is analyzed 
using multiple numerical tools. The chemical equilibrium compositions and applications (CEA) code is 
used to compute the flow field inside the engine nozzle. The plume expansion into ambient atmosphere is 
simulated using an axisymmetric space-time conservation element and solution element (CE/SE) Euler 
code, a computational fluid dynamics (CFD) software. The thermal analysis including both convection 
and radiation heat transfers from the hot gas inside the engine nozzle and gas radiation from the plume is 
performed using Thermal Desktop. Three SM configurations, Lockheed Martin (LM) designed 604, 605, 
and 606 configurations, are considered. Design of multilayer insulation (MLI) for the stowed solar arrays, 
which is subject to plume heating from the main engine, among the passive thermal control system 
(PTCS), are proposed and validated. 


1. Introduction 

The CEV service module will experience heating from the rocket engine when the engine is firing 
during phases of the mission, such as the circularization burn to orbit and the trans-earth injection. The 
main engine and its plume could get extremely hot and will radiate significant heat to its adjacent 
components on the SM. The stowed solar panel is next to the engine nozzle and becomes the major 
concern. Thermal analysis is performed to compute the temperature on the solar panel due to the heat 
from both the main engine and its plume. Proper thermal MLI blanket design for the solar panel is also 
proposed and validated to keep the panel under its maximum safe operating temperature limits. 

In the following, the simulation of the flow field inside the engine nozzle and its plume expansion 
into ambient atmosphere are described first, which is followed by the thermal analysis performed using 
Thermal Desktop. 

2. Simulation of Flow Inside the Engine Nozzle and Its Plume 

As the first step, the CEA code is used to compute the engine performance. The flow properties, such 
as the pressure p , density p, velocity V, temperature T 7 , and mole fractions of the exhaust gas at the 
injector and the end of the combustor, nozzle throat, and nozzle exit, are also obtained. The 
thermodynamic properties of the exhaust gas provided by CEA code are used to compute the heat transfer 
coefficient inside the engine nozzle with the empirical equations in (ref. 1). The flow properties at the 
nozzle exit and the thermodynamic properties are used in the CFD simulation of the plume expansion into 
the ambient air. The gas radiation from the hot exhaust gas is computed using a simplified engineering 
approach described in (ref. 2). Each topic will be described in detail and numerical results are reported as 
follows. 
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2.1 Service Module Main Engine Performance 

The thermodynamic and flow properties of the flow inside the main engine are computed using CEA 
code. Given the key engine parameters, such as the chamber pressure, oxidizer and fuel names and their 
ratio, chamber subtraction and expansion area ratios, propellant mass flow rate as the CEA inputs, the 
CEA computes the properties at different locations inside the thrust chamber and the results are listed in 
table 1 . At the nozzle exit, the mole fractions of the exhaust species based on chemical equilibrium are 
33.56% of H 2 0, 13.6% of C0 2 , 3.22% of CO, 32.8% of N 2 , and 16.9% of H 2 , which will be used in the 
plume radiation calculation. 


TABLE 1.— ENGINE PERFORMANCE (CEA RESULTS) 



Injector 

Combustion 

end 

Throat 

Exit 

p, BAR 

8.6 

8.07 

4.81 

0.00457 

r,K 

3137 

3124 

2967.8 

915.6 

p, kg/m 3 

0.7092 

0.6673 

0.4237 

0.001347 

Sonic velocity, m/s 

1175.6 

1172.7 

1137.3 

651.8 

Mach number 

0 

0.244 

1 

5.082 


2.2 Convective Heat Transfer Inside the Engine Nozzle 

A steady-state one-dimensional model for isentropic flows of an ideal gas through a converging- 
diverging nozzle is used to compute the flow variables of the hot gas inside the engine nozzle along the 
axial direction. The convective heat transfer coefficient that is a function of Reynolds and Prandtl 
numbers based on empirical equations (ref. 1) can be computed. Figure 1 show the computed heat transfer 
coefficient and temperature of hot gas along the flow direction, which will be used later in the thermal 
analysis. 


2.3 Engine Plume Expansion Into the Ambient 

The exhaust plume is simulated using a CFD code, the CE/SE two-dimensional/axisymmetric Euler 
solver. An ideal gas is assumed and the exhaust gas properties at the nozzle exit obtained from CEA is 
used. Two altitudes (100 and 400 km) at the angle of attack a (= 0°) that is the angle between the vehicle 
velocity and thrust vectors are considered. At the inlet of the computational domain (nozzle exit), we 
define 


p = 0.0045 atm, T = 915.6 K, p=1347kg/m 3 , w = 3.312km/s, v = 0.0 km/s 




x(m) x ( m ) 

(a) Static temperature (b) Heat transfer coefficient 

Figure 1 . — Heat transfer coefficient and hot gas temperature inside the nozzle. 
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It should be pointed out that the nozzle cone half angle 20° at the exit of the chamber is not included 
here. The ambient conditions are 


p = p a , p=p a , u = V cosa, v = V sin a 

where p a and p a are the pressure and density of the ambient air, respectively, and V= 7.8 km/s is the 
vehicle velocity. For the 400 km altitude, the CFD code becomes unstable due to numerical issues related 
to very small density and pressure at that altitude. The code can produce converged plume results up to 
100 km altitude. To approximate the ambient air condition at 400 km altitude and above, the pressure and 
density of the ambient air at the altitude of 100 km with zero air velocity are used. It results in zero 
dynamic pressure ( p a V 2 ) of the ambient air and gives a reasonable approximation to the ambient condition 
at 400 km altitude and above since p a V 2 (= 0.00017 N/m 2 ) at the altitude of 400 km is very small. The 
ambient air conditions at different altitudes are listed in table 2. 


TABLE 2.— THE AMBIENT AIR CONDITIONS AT DIFFERENT ALTITUDES 


Altitude, 

Temperature, 

Pressure, 

Density, 

m 

K 

atm 

kg/m 3 

100,000 

195.08 

3.16E-07 

5.60E-07 

200,000 

845.56 

8.36E-10 

2.54E-10 

400,000 

995.83 

1.43E-1 1 

2.80E-12 


The computed numerical results are plotted in figures 2 and 3. In figure 2, logi 0 p, logioP, temperature 
contours and plume relative to CEV are plotted for 100 km altitude at a = 0°. The corresponding plume 
pattern for the altitude of 400 km and above is plotted in figure 3. The air shock, exhaust shock, and 
air/exhaust mixing are seen in the plume at the altitude of 100 km, which agrees with some plume 
patterns shown in (ref. 3). The major mass is within the plume intrinsic core near the nozzle exit. No 
shock waves are observed in the plume and the plume expands more at 400 km altitude. For the vacuum 
case, the exhaust gas expands to reach the Prandtl-Meyer limit that is 103° (ref. 4) including the half cone 
angle of 20°. However, the major mass is still within the intrinsic core. The altitude has negligible effect 
on the intrinsic core. 


2.4 Gas Radiation Computation 

Considering the gas radiation, strong emission by carbon dioxide (C0 2 ) and water vapor (H 2 0) (polar 
molecules) is in infrared (IR) region and weaker emission is in the ultraviolet (UV) and visible ranges. 
Non-polar gases such as 0 2 , N 2 , and H 2 , do not emit radiation and are transparent to incident thermal 
radiation. The approach for computing the gas radiation presented in (ref. 2) is used here and will be 
described briefly in the following. 

The radiation heat flux q = A s s g o , where A s is the surface area, a is Stephan-Boltzmann constant, s g 

is the gas emissivity, correlated in terms of gas temperature T g , total pressure p of the gas, partial pressure 
p g of the radiating species, and mean beam length Le = 4 VJA S with V s being the volume. Furthermore 

Eg = E c + E w — As 

where e w is the emissivity of H 2 0, e c is the emissivity of C0 2 , and As is the correction term for H 2 0 and 
C0 2 mixing. The partial pressure of C0 2 and H 2 0 are computed as 

Pc 02 = P mco2 , Pmo = P m mo 

where m c 02 and m H2 0 are the mole fractions of C0 2 and H 2 0, respectively. The data of s w , s c , and As can 
be found in (ref. 2). 


NAS A/TM— 2007-2 1 5049 


3 




>. 


>. 





(d) Plume relative to CEV 


Figure 2. — Plume at the altitude of 100 km with a = 0°. 
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For the flow inside the engine nozzle, since the flow properties change drastically along the flow 
direction, the hot gas is split into three parts along the axial direction x to calculate the radiation heat flux. 
Each part is defined using its starting and ending x locations, denoted as x^ and x e , respectively. Table 3 
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lists the geometry of the three parts and the corresponding results of the heat flux. Note that the 
regenerative cooling is not included here. 


TABLE 3.— RADIATION HEAT FLUX INSIDE THE NOZZLE 


[Xs, X e ], 

m 

Height, 

m 

Volume, 

m 3 

Surface area, 
m 2 

Le , 
m 

PcoiLe, 

atm-ft 

PmoLe, 

atm-ft 


T* 

K 

kw 

[0, 0.5] 

0.5 

0.1 

0.864 

0.466 

0.454 

1.122 

0.207 

2050 

179 

[0.5, 1.0] 

0.5 

0.417 

1.8059 

0.92456 

0.011 

0.0275 

0.0074 

1050 

0.92 

[1.0, 1.8] 

0.8 

1.7171 

4.52 

1.52 

0.00427 

0.0105 

0.0024 

950 

0.5 


The so-computed radiation heat flux for the plume is listed in tables 4 and 5 for 100 and 400 km 
altitudes, respectively. The volume-averaged gas pressure and temperature are used in the calculation. 
Since the regenerative cooling in the engine nozzle is not considered, the actual exhaust gas temperature 
will be lower. For the 100 km altitude, the plume core is split into three right circular cones (shown in 
fig. 2(b)) based on the pressure contour. 


TABLE 4.— MAIN ENGINE EXHAUST PLUME RADIATION HEAT FLUX (100 km, a - 0°) 


[ri,r 2 \, 

m 

Height, 

m 

Volume, 

m 3 

Surface area, 
m 2 

Le , 
m 

PcoiLe, 

atm-ft 

PmoLe, 

atm-ft 


T 

K 

i b 
kw 

[LI, 5.6] 

10.1 

410 

233 

4.89 

0.0011 

0.00275 

0.0008 

850 

5.5 

[5.6, 7.0] 

5.4 

676 

221 

5.7 

0.000688 

0.0017 

0.0005 

600 

0.8 

[7.0, 8.5] 

7 

1325 

348.6 

7.265 

0.00029 

0.0007 

0.00025 

500 

0.3 


TABLE 5.— MAIN ENGINE EXHAUST PLUME RADIATION HEAT FLUX (400 km ALTITUDE AND ABOVE) 


[ri,r 2 ], 

Height, 

Volume, 

Surface area, 

Le , 

PcoiLe, 

Pmo^e, 

E g 

T 

<L 

m 

m 

m 3 

m 2 

m 

atm-ft 

atm-ft 

K 

kw 

r4,ioi 

10.8 

1764.3 

543.4 

7.774 

0.00066 

0.001635 

0.0004 

850 

6.43 

[10,15] 

13.2 

6566 

1108.6 

12.333 

0.000456 

0.001128 

0.00018 

550 

1.03 


Each right circular cone is defined using its top and bottom radius denoted as r 7 and r 2 and the height. 
Thus the corresponding surface area and volume can be computed. The radiation due to the interaction 
between the exhaust and ambient air is not calculated here. For 400 km and above altitudes, the plume 
core is split into two right circular cylinders (shown in fig. 3(b)). The radiation due to the interaction 
between ambient air and exhaust gas is negligible due to the high altitude. 

The radiation from both the rocket nozzle and the expanding exhaust hot gas are being considered in 
the thermal analysis. The gas radiation is a function of the combustion products and their density and 
pressure. As a result, over 95% of the heat produced in the exhaust plume is within the intrinsic core of 
the exhaust plume. The gas that expands toward the Prandtl-Meyer expansion limit has very low pressure 
and density and does not generate significant heat. However, the radiation from the hot engine nozzle 
produces major heating. As a result, the radiation shielding for the adjacent components on the SM might 
be required based on further analysis. 

3. Thermal Analysis on the Service Module Components 

The thermal analysis is performed using Thermal Desktop. The major components in the model are 
the SM main engine, solar panels in a stowed position, and the radiator panel. Inside the main engine, 
both convection and radiation boundary conditions are imposed. The exhaust plume was modeled as a 
series of connected parabolas that are defined based on the radius and height listed in Table 3, and the 
corresponding heat flux is applied to each section of the parabola. The total heat radiated from the 
exhaust plume was 6.6 kW of heat for the 100 km altitude and 7.46 kW for the 400 km altitude. Thermal 
Desktop was used to calculate the view factors between the geometry and the resulting temperatures. The 
stowed solar panel whose optical properties used here is an emissivity of 0.9 and solar absorptivity of 0.9. 
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The ten layers of MLI insulation comprises a 3.0-mil-thick black Kapton outer layer, and four layers of 
double-aluminized Kapton, and five layers of double-aluminized Mylar with Dacron netting spacers. Its 
effective emmisivity is 0.05 and solar absorbtivity is 0.14 (ref. 5). The five layer design comprises an 
outer layer of double-aluminized Kapton and four layers of double-aluminized Mylar with Dacron 
spacers. Its effective emmisivity is 0.1 and solar absorbtivity is 0.14 (ref. 5). Three configurations of CEV 
SM (LM 604, 605, and 606) are considered. For each configuration, three cases including no insulation, 5 
layers of MLI, and 10 layers of MLI are studied by changing the optical properties on the solar panel with 
its corresponding insulation. On each surface of the solar panel, two mesh elements are used in the model 
for the 604 configuration; three mesh elements are used in the model for the 605 and 606 configurations. 
More mesh elements can be used if it is necessary. The ambient temperature is assumed to be at 0 K. 

For the CEV 604 configuration, the computed temperature contour at t = 0.5 hr with five layers of 
MLI is plotted in figure 4 to show a representative temperature distribution on the thermal model. 
Furthermore, the time history of the temperature at the center of two mesh elements on the solar panel 
surface is plotted in figure 5 for different thermal insulations. It can be seen that, with no insulation, it 
takes about 0.2 hr for the solar panel to reach the steady state at 650 °F, which exceeds the safe-operating 
limit of the solar panel. The temperature reaches 575 °F for using a 5 -layer MLI blanket and 420 °F for 
using a 10-layer MLI blanket at t = 0.5 hr. The steady state has not been reached yet. 

For the CEV 605 configuration, the computed temperature contour at t = 0.5 hr with 10 layers of MLI 
is plotted in figure 6, while the time history at the center of the three mesh elements on the solar panel 
surface with different insulations is plotted in figure 7. Similarly, the maximum temperature on the solar 
panel could reach 680 °F within 0.3 hr if there is no insulation on the solar panel. The maximum 
temperature could drop to 625 °F using 5 layers of MLI and 425 °F using a 10-layer MLI blanket at t = 
0.5 hr. 

For the CEV 606 configuration, the computed temperature contour at t = 0.5 hr with five-layer MLI is 
plotted in figure 8 and the corresponding time history is plotted in figure 9. The solar panel is farther 
away from the engine nozzle and most of the area is blocked by the radiator panel compared with the 604 
and 605 configurations. Without a MLI blanket, the maximum temperature is 375 °F at t = 0.3 hr and 
reaches steady state. With a 5-layer MLI, the maximum temperature is 230 and 145 °F for 10 layers of 
MLI at t = 0.5 hr. 



Tenpercrture CFJ. Tine = 0, 5 hr 


case3£, sav 


Figure 4. — Temperature contour at t = 0.5 hr for the 604 
configuration with 5-layer MLI blanket. 
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Time, hr 



Time, hr 



Time, hr 

Figure 5. — Time history of temperature at the center of two 
mesh elements on solar panels for the 604 configuration, 
(a) With no insulation, (b) Five-layer MLI, (c) Ten-layer MLI. 



c nse3g , sav 


Figure 6. — Temperature contour at t = 0.5 hr for the 
605 configuration with 10-layer MLI blanket. 
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Time, hr 



Time, hr 



Time, hr 

Figure 7. — Time history of temperature at the center of 
three mesh elements on solar panels for the 605 
configuration, (a) With no insulation, (b) Five-layer MLI, 
(c) Ten-layer MLI. 



Figure 8. — Temperature contours at t = 0.5 hr for 
the 606 configuration with 5-layer MLI blanket. 


NAS A/TM— 2007-2 1 5049 


9 




Time, hr 



Time, hr 



Time, hr 

Figure 9. — Time history of temperature at the center of 
three mesh elements on solar panels for the 606 
configuration, (a) With no insulation, (b) Five-layer MLI, 

(c) Ten-layer MLI. 

In summary, the CEV 606 configuration provides the best thermal environment from the main engine 
for the solar panel. Both the 604 and 605 configurations have similar thermal conditions and need more 
layers of MLI blanket. 


References 

1. F.P. Incropera and D.P. DeWitt, “Fundamentals of heat and mass transfer,” 5 th Ed., John Wiley & Sons 
Inc., 2002. 

2. R. Siegel and J.R. Howell, “Thermal radiation heat transfer,” 3 rd Ed., Hemisphere Publishing 
Corporation, Washington, 1992. 

3. F.S. Simmons, “Rocket exhaust plume phenomenology,” p. 25, Aerospace Corporation, El Segundo, 
CA, 2000. 

4. M.A. Saad, “Compressible fluid flow,” Prentice-Hall Inc., 1995. 

5. D.G. Gilmore, “Spacecraft thermal control handbook,” Volume I: fundamental technologies, 2 nd Ed., 
The Aerospace Press, El Segundo, CA, 2002. 


NAS A/TM— 2007-2 1 5049 


10 


REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the 
data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this 
burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. 
Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB 
control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-YYYY) 
01-11-2007 

2. REPORT TYPE 

Technical Memorandum 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

Thermal Analysis on Plume Heating of the Main Engine on the Crew Exploration Vehicle 
Service Module 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Wang, Xiao- Yen, J.; Yuko, James, R. 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

WBS 931164.03.01.02.01 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E- 16260 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSORING/MONITORS 
ACRONYM(S) 

NASA 

11. SPONSORING/MONITORING 
REPORT NUMBER 

N AS A/TM-2007 -2 15049 

12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category: 34 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390 

13. SUPPLEMENTARY NOTES 

14. ABSTRACT 

The crew exploration vehicle (CEV) service module (SM) main engine plume heating is analyzed using multiple numerical tools. The 
chemical equilibrium compositions and applications (CEA) code is used to compute the flow field inside the engine nozzle. The plume 
expansion into ambient atmosphere is simulated using an axisymmetric space-time conservation element and solution element (CE/SE) 
Euler code, a computational fluid dynamics (CFD) software. The thermal analysis including both convection and radiation heat transfers 
from the hot gas inside the engine nozzle and gas radiation from the plume is performed using Thermal Desktop. Three SM configurations, 
Lockheed Martin (LM) designed 604, 605, and 606 configurations, are considered. Design of multilayer insulation (MLI) for the stowed 
solar arrays, which is subject to plume heating from the main engine, among the passive thermal control system (PTCS), are proposed and 
validated. 

15. SUBJECT TERMS 

Thermal analysis; Plume heating 

16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

UU 

18. NUMBER 
OF 

PAGES 

16 

19a. NAME OF RESPONSIBLE PERSON 

STI Help Desk (email:help@ sti.nasa.gov) 

a. REPORT 

U 

b. ABSTRACT 

U 

c. THIS 
PAGE 

U 

19b. TELEPHONE NUMBER (include area code) 

301-621-0390 


Standard Form 298 (Rev. 8-98) 
Prescribed by ANSI Std. Z39-18 















